Resolving the plasma profile via differential single inclusive suppression 
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The ability of experimental signatures to resolve the spatio-temporal profile of an expanding quark 
gluon plasma is studied. In particular, the single inclusive suppression of high momentum hadrons 
versus the centrality of a heavy-ion collision and with respect to the reaction plane in non-central 
collisions is critically examined. Calculations are performed in the higher twist formalism for the 
modification of the fragmentation functions. Radically different nuclear geometries are used. The 
influence of different initial gluon distributions as well as different temporal evolution scenarios on 
the single inclusive suppression of high momentum pions are outlined. It is demonstrated that the 
modification versus the reaction plane is quite sensitive to the initial spatial density. Such sensitivity 
remains even in the presence of a strong elliptic flow. 



PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 



The goal of ultra-high-energy heavy-ion collisions is 
the creation and study of strongly interacting matter, 
heated past a temperature beyond which confinement 
can no longer be expected Q. The Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory 
(BNL), has provided a wealth of data which call for de- 
tailed theoretical modeling and study of the produced 
matter At present, computations of the wide vari- 
ety of bulk obscrvablcs 0. |j. have led to a picture 
that the Quark-Gluon Plasma (QGP) formed at RHIC 
assumes the dynamical properties of an ideal fluid 0. 
The study of hard probes complements this view by 
ascribing considerable opaqueness to the produced mat- 
ter towards the passage of hard jets through it. This set 
of observables, when taken in its entirety, has led to the 
suggestion that the matter formed at RHIC is strongly 
interacting @. 

The object of this letter is to study the feasibility of 
correlating the properties of the plasma derived from in- 
vestigations in the bulk sector with those gleaned from 
the hard probe sector. Full 3+1 dimensional hydrody- 
namic simulations have performed very well in com- 
parison to experimental data in the soft sector. In 
essence, these simulations make predictions regarding the 
space-time profile (STP) of the expanding matter, but re- 
quire, as input, an initial spatial density profile and are 
sensitive to the ansatz used. The possibility of testing 
such a space-time profile via an independent set of mea- 
surements afforded by the modification of the fragmen- 
tation of hard jets as they pass through the dense matter 
forms the focus of the current study. Such a study is 
similar, in spirit, to Ref. 0] which explored the effect 
of the longitudinal flow on the modification of hard jets. 
This complementary effort will assume boost invariancc 
and will be restricted to mid-rapidity. 

There exists a wide variety of observables involving 
particles at high transverse momentum px which may 
offer insight regarding the space-time profile of the mat- 
ter produced in a heavy-ion collision: single inclusive ob- 
servables (both differential and integrated) 0, 0] , dou- 
ble inclusive hadronic observables (both near and away 



side) mm> 

as well as photon hadron correlations [Tlj . 
In this letter, the ability of single inclusive observables to 
reveal the space-time profile of the plasma will be criti- 
cally examined. It has already been pointed out that the 
nuclear modification factor Raa as a function of px for 
the most central event is not very sensitive to the detailed 
nature of the initial gluon density or the time evolution 
of this density 01- The induction of such sensitivity 
requires the use of more differential probes. In this let- 
ter, the nuclear modification versus the reaction plane 
in non-central collisions will be examined. A simplified 
methodology will be followed: While both the production 
cross section of a jet (which depends on the number of bi- 
nary collisions A^in) as we H as the density of the plasma 
at a certain transverse location (which is more closely re- 
lated to the number of participants A part ) depend on the 
nuclear density profile, these will be disassociated from 
one another. The nuclear density profile relevant to the 
estimation of the number of initial binary collisions may 
be determined independently in pA collisions. For the 
current treatment, this will be approximated by a hard 
sphere distribution and left unchanged. As the focus is 
on the dependence of the jet modification on the profile 
of the produced matter in the final state, different nu- 
clear densities and temporal dependences will be used as 
inputs to produce a space-time modulation of the pro- 
duced matter. In this way, a direct measure of the den- 
sity profile-dependent modification on an identical set of 
produced jets will be explored. 

The phenomenology of particle production at very high 
energies in pp collisions is greatly simplified by the fac- 
torization theorems of QCD m> These predict that at 
high transverse momentum px (and as a result large vir- 
tually Q 2 ) the single inclusive cross-section achieves a 
factorized form. For high enough transverse momentum 
Pt > 7GeV, further simplifications arise in the case of 
nuclear collisions, i.e., one may ignore a variety of initial 
and final state nuclear effects, such as recombination, in- 
trinsic transverse momentum and the Cronin effect. At 
such high energies these are known to be small com- 
pared to the dominant effect of medium induced energy 
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loss which leads to the modification of the fragmentation 
function [TtI ]. Such effects may have an interesting den- 
sity dependence and will reappear in any effort to extend 
the effects discussed herein to lower pr- Assuming a fac- 
torization of initial and final state effects, the differential 
cross-section for the production of a high pt hadron at 
midrapidity from the impact of two nuclei A and B at 
an impact parameter between 6 m i„, b max is given as, 



da 



AB 



dyd 2 p T 



K d 2 b d 2 rt A (r + b/2)t B (r-b/2) 



x Jdx a dx b G*(x a ,Q 2 )Gi(x bl Q 2 ) 

y derated D*(z,Q 2 ) 
di 7TZ 



(1) 



AD(z,Q 2 ,r) = g/^/ ^P^M^sCa 



x T M {b,f,x a ,x b ,y,l ± )D i - Q 



l 2 ± N c t A (f+b/2)t B (r-b/2) 



x Gi(x a )G*(x b )§ 
dt 



+ v.c. 



(3) 



In the above equation, l± is the transverse momentum 
of the radiated gluon (quark) which leaves a momen- 
tum fraction y in the quark (gluon) denoted as parton i 
(P q ^i(y) is the splitting function for this process) which 
then fragments leading to the detected hadron. The v.c. 
refers to virtual corrections. The factor T M originates 
from the higher twist matrix element, which encodes the 
information of rescattering off the soft gluon fields, 



where, G^(x a ,Q 2 )Gf : (x b ,Q 2 ) represent the nuclear par- 
ton distribution functions. These aregiven in terms of 
the shadowing functions S A (x a , Q 2 ) [l8| and the parton 
distribution functions in a nucleon G a (x a ,Q 2 ) as 



T = t A (f+b/2)t B (f-b/2)Gi(x a )Gf(x b ) 



G a (x a , 



S A (x a ,Q 2 )G a (x a ,Q 2 ). 



(2) 



In Eq. t A (r),tB(r) are the thickness functions of nu- 
clei A and B at the transverse location r = {x,y). As 
mentioned previously, throughout this letter, such thick- 
ness functions will be evaluated using the simple and 
analytically tractable hard sphere densities (p(x,y,z) = 
p$0(B 2 A — x 2 — y 2 — z 2 )). In this formulation, s, t, u re- 
fer to the Mandelstam variables of the internal partonic 
process. Unless, stated otherwise, the variables without 
the hats refer to the variables of the full process. The 
two participating partons in the initial state are referred 
to as a, b while the final state partons are referred to as 
c, d. The factor K ~ 2 accounts for higher order contri- 
butions. 

The most important element in Eq. is the medium 
modified fragmentation function D(z,Q 2 ) expressed as 
the sum of the leading twist vacuum fragmentation func- 
tion and a correction brought about by rescattering of the 
struck quark in the medium i.e., D = D + AD The 
vacuum fragmentation functions D(z, Q 2 ) are taken from 
Ref. [2(|. In the collinear limit, the modification is com- 
puted by isolating corrections, suppressed by powers of 
Q 2 , which are enhanced by the length of the medium |2lj. 
At next-to- leading twist, the correction for the fragmen- 
tation of a quark has the expression (generalized from 
deep-inelastic scattering (DIS) [T^.[22T| ). 



da 
dt 



x / d(x g p g {x gi n( + r)(2-2cos(r) L ()). (4) 
Jo 

In the above equation, the factor tjl = l\l\2p>Ty{X ~ V)\- 
Where, px represents the transverse momentum of the 
produced parent jet. It should be pointed out that in 
the higher twist formalism used here, no assumption is 
made regarding the prevalent degrees of freedom of the 
produced matter. 

Breaking with usual practice, the jet direction (n) is 
chosen as the x-axis in the transverse (x, y) plane (the 
z-direction is set by the beam line); the angle of the re- 
action plane vector b is measured with respect to this 
direction. The distance travelled by the jet in this di- 
rection prior to scattering off a gluon is denoted as £. 
The gluon's forward momentum fraction in denoted as 
x g . One assumes, following Ref. j2^| that the relevant 
values of small enough that x g p g (x g ,n() is almost 

independent of x g and the combination will henceforth 
be referred to as the density p. This represents the one 
unknown in the above set of equations and in the higher 
twist formalism. It is related to the more familiar quan- 
tity q by the simple equation |23l |. 



q(f+hC) 



4Tr 2 a s C F 
N3 - I 



p(x g ,r + h(). 



(5) 



This is the q at a location (and time) r+nC in the plasma, 
and not the averaged value. 

Ignoring the modification in the fragmentation func- 
tion {AD) and the factors of nuclear shadowing in 
Eq. produce the result of binary scaling the p — p 
scattering cross section and produce the denominator of 
the nuclear modification factor, 
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Raa = 



TAA{b 



, \ d<r™{p T ,y) 
mm ) "max ) dyd 2 p T 



(0) 



The phenomenological input that is required to un- 
derstand the variation of the nuclear modification factor 
with centrality or pt is the space time dependence of 
the gluon density p(x,y,z,r). In the remainder, the fo- 
cus will be on observables at midrapidity and thus the 
z-coordinatc will be ignored; different forms of the gluon 
density profile will be invoked and their effect on the 
nuclear modification factor versus centrality and the re- 
action plane will be elucidated. 

The initial spatial profile of the plasma density depends 
on the nuclcon density brought in by the colliding nuclei. 
It, most likely, should be a smooth, monotonously drop- 
ping function of the transverse coordinates (x,y), the 
time r as well as the impact parameter of the collision 6. 
The function must approach zero as any of these parame- 
ters are increased from their minimum values. A general 
functional form which satisfies the above constraints is 



p(r, t; b) = p Q 



to N wn {t f A (r + b/2),t f A (f- 6/2)) 



N wn (t f A (o),t A m 



(7) 



where, N wn is the wounded nucleon participant density 
HQ, and t f A (f±b/2) is an input thickness function that 

may be varied to obtain different spatial profiles p(r, b, t). 
The superscript / on the thickness functions is meant 
to differentiate it from the thickness functions used to 
generate the number of binary collisions in the initial 
state. The time dependence of the density for the first 
set of computations will be assumed to be generated by 
a pure Bjorken expansion (which leads to a 1/t depen- 
dence in Eq. J7J). It may be immediately verified that 
even within the apparently constraining form of Eq. iJJJ, 
radically different space-time profiles may be obtained by 
using different input thickness functions. The profiles are 
normalized by the maximum density achieved at r = 
in the most central collision at b = (the denominator of 
Eq. 10). In such a form, the initial spatial modulation 
is contained entirely in the normalized modulation factor 
M(x,y) i.e., p(x,y) = po m M(x,y). By definition, the 
maximum value of M(x,y) = 1. The unknown factor, 
PqTq in the above equation is set in all cases by assuming 
that the Raa for the most central collisions at a pion 
Pt = 8 GeV be 0.2. Once so normalized, the variation of 
the Raa with pr for all the different profiles considered 
for the most central collision is remarkably similar. The 
above model for the density may then be used to predict 
the Raa at a lesser centrality and also the variation of 
the Raa with the reaction plane. 

In an effort to demonstrate the sensitivity of the nu- 
clear modification factor to the space-time profile, two 
rather disparate input thickness functions t A are used: 



the hard sphere and Gaussian. The hard sphere thick- 
ness function is given as, 



t f A (r, 6/2) = C m ^/R A -(x± b x /2Y - (y ± by/2)*. (8) 

The focus in this letter will be on the peripheral events 
with a centrality between 40% — 60%, there exists a con- 
siderable azimuthal asymmetry in the gluon distribution 
in such an event for a hard sphere distribution. The re- 
sults of such a calculation are shown in Fig. ^ I n this 
figure, the solid black line represents the Raa integrated 
over all emission angles. The other lines represent the dif- 
ferential Raa as a function of the angle of the jet with re- 
spect to the reaction plane (cf>) i.e., Raa{^4>)- Due to the 
large azimuthal asymmetry in the initial gluon density, 
leading to jets facing different amounts of matter depend- 
ing on emission angle, there exists a large asymmetry as 
well in the R AA {&4>). The plain dashed line (0 = 0-15°) 
represents emission in the reaction plane. These repre- 
sent jets which pass through the least amount of matter 
and as a result exhibit the lowest amount of modification. 
The dashed line with triangles (<j> = 75 — 90°) represents 
emission perpendicular to the plane. These pass through 
the maximum amount of matter and thus show the most 
modification. The lines in between represent emission at 
intermediate angles. This figure clearly demonstrates the 
sensitivity of the nuclear modification factor versus the 
reaction plane to the azimuthal asymmetry of the gluon 
distribution. 
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FIG. 1: (Color online) Raa at 40-60% centrality (solid black 
line) and the variation of Raa as a function of the angle cj> 
with respect to the reaction plane, assuming the initial gluon 
density generated by the overlap of two hard-sphere geome- 
tries. The evolution of the density is assumed to be purely 
due to Bjorken expansion. 

To put the above statement on a firm footing, the sen- 
sitivity of the Raa{&4>) is tested against a gluon dis- 
tribution which docs not posses as large an azimuthal 
asymmetry as in the case above. This is achieved by 
using a Gaussian density as the input thickness function, 
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dependence of r and e may be obtained in a toy model: 



t{(f,b)=C G e 



\f±b/2[ 



(9) 



where Rq is an appropriately chosen Gaussian radius. As 
would be expected, such an initial gluon density leads to 
a very different variation of the Raa with respect to the 
reaction plane. This is plotted in Fig.El The variation of 
the nuclear modification factor with <f> is much reduced as 
compared to that from a hard sphere thickness function. 
This difference testifies to the effectiveness of observables 
such as the Raa(^4>) as probes of the initial gluon den- 
sity profile. For the initial number of binary collisions, 
the use of a hard sphere density is continued. This en- 
sures that all jets are produced in the same, finite, well 
defined, volume. 
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FIG. 2: (Color online) Same as Fig. with the hard sphere 
density replaced with a Gaussian profile. 

The variation of the Raa with <f> may also be sensi- 
tive to the temporal evolution of the plasma. In both 
of the preceding calculations, the temporal dependence 
was taken to be that due to Bjorken expansion. In real 
heavy-ion collisions, there is a considerable amount of 
transverse (radial and elliptic) flow of the matter caused 
by the large pressures built up in the dense matter. To 
have a completely realistic description of such a scenario 
requires that the density modulation factor M(x, y, t) 
be extracted from full three dimensional hydrodynami- 
cal simulations [24|. In this effort, a simple scaling form 
will be employed: the two transverse coordinates are re- 
placed by x/[r(r)e(T)] and ye(r)/r(r), where r(r) de- 
scribes the radial expansion and e(r) generates the az- 
imuthally asymmetric expansion. It should be noted that 
there already exists an asymmetry in the density dis- 
tribution from Eq. J7J) where the hard sphere thickness 
function (Eq. |SJ is used as input to calculate the gluon 
density. The factors r, e are meant to introduce a further 
time dependent modulation. Thus at r = 0, both these 
factors must be set to unity. A simple form for the time 



r(r) = H t and e(r) 2 = 1 

Ra 



2Ra 



t. (10) 
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FIG. 3: (Color online) Same as Fig. Q including transverse 
flow. 

In the above, the radial flow velocity v is taken to be 0.4 
from Ref. The impact parameter dependent asym- 
metric velocity v e is taken to be such as to restore the 
azimuthal symmetry to the boundary of the expanding 
plasma by r = 2R A i.e., v e = y ma x{b) / x max (b) - 1. The 
maximum x and y coordinate of the overlap at r = is 
represented as x max , y max (note that the use of the hard 
sphere thickness function as input allows a clear identifi- 
cation of an x m ax and y max )- The integral of the density 
over the three volume is kept fixed with the replacement 



r(r) 2 ' 



(11) 



The result of such a model calculation is shown in Fig. [3] 
As is clearly demonstrated, the large asymmetric flow, 
leads to an Raa versus reaction plane with a lesser spread 
as a function of <f> than in Fig. ^ The ordering of the 
contributions remains unchanged. Thus even with a re- 
alistic amount of radial and elliptic flow, the Raa(A4>) 
still maintains a sensitivity to the initial asymmetry. If 
the asymmetric expansion is increased, this leads to a 
further reduction in the spread. This is shown in Fig. 0] 
where the Raa, along with the Raa in reaction plane 
and out of plane are plotted as a function of f e , an over- 
all multiplicative factor used to modify the asymmetric 
velocity v e . In the figure, f e = corresponds to the case 
of Fig. with a radial flow, i.e., with no elliptic flow. The 
case of v e = 1 corresponds to the case of Fig- El *- e -> with 
elliptic flow parametrized from Ref. 0. Higher values 
of / e represent multiples of the asymmetric velocity used 
in Fig. |21 As expected, raising the asymmetric velocity 
leads to a decrease in the spread of the Raa as a function 
of the reaction plane. It should be pointed out that in 
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both Figs. ^ an d El the final state gluon density profile 
starts out with a considerable azimuthal asymmetry as a 
hard sphere thickness function is used as input in Eq. JjJ. 
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FIG. 4: (Color online) Effect of asymmetric flow on the spread 
of Raa as a function of the reaction plane at pion pr = 10 
GeV. Fhe solid line is the angle integrated Raa the dashed 
lines with (without) triangles is the Raa with the jet emission 
angle out of (in) the reaction plane. 

The calculations presented above offer sufficient evi- 
dence that differential probes such as the Raa versus 
the reaction plane offer deeper insight into the space- 



time profile of the matter produced at RHIC. Sharper 
eccentricities lead to large variations of the Raa with 
respect to the emission angle. The effect of large ellip- 
tic flow is discernable. It is true that such observations 
alone will not distinguish the effect of a large elliptic flow 
from that of a lower eccentricity in the initial condition. 
To discern between such facets will require multiple dif- 
ferential probes to be compared in unison. The nuclear 
modification factor versus the reaction plane carries more 
information than the azimuthal asymmetry of high px 
hadrons (the v-i at high pt)- The reader will note that 
the V2 obtained from the Raa(^) (Figs. 1113)1 is insensi- 
tive to the overall magnitude of the Raa- 

The current exercise is meant to serve as a first attempt 
in the resolution of the spatio-temporal profile of the hot 
plasma. As such, a variety of approximations were made, 
kr broadening was ignored, and the number of initial 
binary collisions was taken from a hard sphere nuclear 
density distribution. This was done in the interest of 
focussing the study on the effect of the plasma profile 
on the variation of Raa with respect to the jet emission 
angle. The inclusion of such effects will result is slender 
modifications to the results outlined in this letter and are 
left for a future, more rigorous comparison with data. 
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